BOUCHAUD WALKS WITH VARIABLE DRIFT 



MANUEL CABEZAS PARRA 



Abstract. We study the possible scaling limits of a sequence of Bouchaud trap models on Z 
with a drift which decays to as we rescale the walks: as the time parameter is rescaled by n, 
the drift decays as n~ a , for some fixed a > 0. Depending on the speed of the decay of the drift 
we obtain three different scaling limits. If the drift decays slowly as we rescale the walks (small 
a), we obtain the inverse of an o-stable subordinator as scaling limit. If the drift decays quickly 
as we rescale the walks (big a), we obtain the F.I.N, diffusion as scaling limit. There is a critical 
speed of decay separating these two main regimes, where a new process appears as scaling limit. 
This new process is a drifted Brownian motion with a random, purely atomic speed measure. 
The critical speed of decay a c of the drift is related to the index a of the inhomogeneity of the 
environment by the equation a c = a/(a + 1). 
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1. Introduction 

The Bouchaud trap model (B.T.M.) is a continuous time random walk X on a graph Q with 
random jump rates. To each vertex x of Q we assign a positive number t x where {t x ) x ^q is an 
i.i.d. sequence such that 

lim u a P[T x > u] = 1 (1) 

U— ¥00 

with a G (0, 1). This means that the distribution of t x has heavy tails. Each visit of X to x G Q 
lasts an exponentially distributed time with mean t x . Let S{k) be the time of the fc-th jump of X. 
(S(k), k G N) is called the clock process of X. Let := X(S(k)) be the position of X after the 
fc-th jump. (Yfc : k G N) is called the embedded discrete time random walk associated to X. 

This model was introduced by J. -P. Bouchaud in 9 and has been studied by physicists as a 
toy model for the analysis of the dynamics of some complex systems such as spin-glasses. More 
precisely, each vertex x of Q corresponds to a metastable state of the complex system, and X 
represents the trajectory of the system over its phase space. One of the phenomena that this 
model has helped to understand is that of aging, a characteristic feature of the slow dynamics of 
many metastable systems. For an account of the physical literature on the B.T.M. we refer to [10] . 

The model has also been studied by mathematicians on different graphs, exhibiting a variety 
of behaviors. In 11 , Fontes, Isopi and Newman analyze the one-dimensional case (Q = Z)and 
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where the walk X is symmetric. They obtain a scaling limit for X which is called the Fontes- 
Isopi- Newman (F.I.N.) singular diffusion. This diffusion is a speed measure change of a Brownian 
motion by a random, purely atomic measure p, where p is the Stieltjes measure associated to 
an a-stable subordinator. Different aging regimes for the one-dimensional case where found by 
Ben-Arous and Cerny in [J]. In higher dimensions (Q = 1 d ,d > 2), the symmetric model has a 
behavior completely different to the one-dimensional shown by Ben Arous and Cerny in 

[5J, and by Ben Arous, Cerny and Mountford in JS]. In these papers, a scaling limit and aging 
results where obtained for X. The scaling limit is called fractional kinetic process (F.K.P) which is 
a time-change of a d-dimensional Brownian motion by the inverse of an a-stable subordinator. In 
PQ and [5] Ben Arous, Bovier and Gayrard obtained aging properties of the model on the complete 
graph. A study of this walk for a wider class of graphs can be found on [7] . For a general account 
on the mathematical study of the model, we refer to [5]. 

The difference between the one dimensional case and the model in higher dimensions can be 
understood as follows. We can express the clock process S(k) of X as S(k) — J2i=o T^e*) where 
the are standard i.i.d. exponential random variables. Thus, the increments of S(k) are the 
depths of the traps (t x ) x& q as sampled by Yfc. In the model in dimensions higher than two, the 
embedded discrete time random walk Yj. is transient (the case d = 2 is more delicate). Thus 
will sample each trap t x a finite number of times. That implies that S(k) does not have long range 
interactions with its past and its scaling limit will be a Markovian process, which is an a-stable 
subordinator. On the other hand, in the one-dimensional symmetric B.T.M., we have that the 
embedded discrete time random walk Yfc is recurrent. Thus will sample each trap t x an infinite 
number of times. In this case, S(k) has long range interactions with its past and its scaling limit 
will be non-Markovian. Furthermore, the clock process S(k) will converge to the local time of a 
Brownian motion integrated against the random measure p. Here p plays the role of a scaling limit 
for the environment (t x ) x q%. 

It is natural to ask if we can find intermediate behaviors between the transient case (d > 1) 
and the recurrent case (d = 1): if we introduce a drift to the one-dimensional B.T.M., note that 
the embedded discrete random walk becomes transient. Thus, intermediate behaviors between the 
transient and the recurrent case might appear when one analyzes a sequence of one-dimensional 
B.T.M.'s with a drift that decreases to as we rescale the walks. In this paper we study this 
question, showing that the speed of decay of the drift sets the long-term behavior of the model and 
exhibiting a sharp phase transition in terms of the type of limiting processes obtained. We next 
describe with more precision the way in which we define the B.T.M. with drift and the results that 
are obtained in this paper. 



BOUCHAUD WALKS WITH VARIABLE DRIFT 



3 



For each e > 0, denote by X e the B.T.M. on Z where the transition probabilities of the embedded 
discrete time random walk are to the right and =^ to the left. We will call this process 
the B.T.M. with drift e. For a > 0, consider a rescaled sequence of B.T.M's with drift n~ a , 
(h a (n)X n (tn))t > 0), indexed by n, where h a (n) is an appropriate space scaling depending on 
a. We will see that as the drift decays slowly (small a), the sequence of walks converges to the 
inverse of an a-stable subordinator, whereas if the drift decays fast (large a) the limiting process 
is the F.I.N, diffusion. As these two posibilities are qualitatively different, we are led to think that 
there is either, a gradual interpolation between these two behaviors as the speed of decay changes, 
or a sharp transition between them as the speed of decay changes. We establish that there is a 
sharp transition between the two scaling limits, that there is a critical speed of decay where a 
new, previously, process appears and that the transition happens at a — a /{a + 1). As the main 
theorem of this paper, we prove that, depending on the value of a, there are three different scaling 
limits: 

• Supercritical case (a < a/(a + 1)). The sequence of walks converges to the inverse of 
an a-stable subordinator. 

• Critical case (a = a /(a + 1)). The sequence of walks converges to a process which is a 
speed measure change of a Brownian motion with drift that we will call the drifted F.I.N, 
diffusion. 

• Subcritical case (a > a/(a+l)). The sequence of walks converges to the F.I.N, diffusion. 

The case a = (contained in the supercritical case), which corresponds to a constant drift, was 
already addressed by Zindy in [15] . 

Let us now make a few remarks concerning the proof of our main theorem. The strategy of 
the proof for the supercritical case is a generalization of the method used in |15j and relies on the 
analysis of the sequence of processes of first hitting times (H^(x);x £ [0,nS]) (S is fixed, b > 0) 
defined as 

H£{x) :=inf{f : X n ~ b (t) >x}. (2) 

We show that these processes (properly rescaled) converge to an a-stable subordinator. From that, 
it follows that the maximum of the walks converges to the inverse of an a-stable subordinator. 
This part of the proof requires some care, because, as we are working with a sequence of walks 
with variable drift, we cannot apply directly the methods used in |15j . It turns out that we have to 
choose b properly to obtain a sequence of walks with the desired drift as we invert the hitting time 
processes. Then, it is easy to pass from the maximum of the walk to the walk itself. In [IT] The 
proof corresponding to the critical case follows the arguments used by There they express 
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rescaled, symmetric one-dimensional B.T.M.'s as speed measure changes of a Brownian motion 
trough a random speed measure. But here we are working with asymmetric walks, so we cannot 
work with the expression used there. To treat the asymmetry of the walks, we use a Brownian 
motion with drift instead of a Brownian motion. That is, we express each walk X n " + ' as a 
speed measure change of a Brownian motion with drift, and then prove convergence of the sequence 
of speed measures to p. The latter is achieved by means of a coupling of the environments. In the 
subcritical case, although we obtain the same scaling limit as in [TT] (a F.I.N, diffusion), again, 
because of the asymmetry of the model, we cannot work with the expression used there. We deal 
with this obstacle using, besides a random speed measure, a scaling function. That is, we express 
the rescaled walks as time-scale changes of a Brownian motion. Then we prove that the scale 
change can be neglected and show convergence of the sequence of speed measures to the random 
measure p. 

The organization of the paper is as follows. In section [2] we give the definition of the model and 
state our main results. There we also give simple heuristic arguments to understand the transition 
at a = a/(a + 1). In section[3]we obtain the behavior for the supercritical case, and in section^ 
we obtain the scaling limit for the critical case. The behavior for the subcritical case is obtained 
in section [5j 

Finally, we would like to mention that while preparing the final version of this article we have 
learned that Theorem (|2.1I) has been independently obtained by Gantert, Morters and Wachtel 
[12] . There, they also obtain aging results for the B.T.M. with vanishing drift. 

Acknowledgements: The author was supported by a fellowship of the National Commission 
on Science and Technology of Chile (Conicyt)#29100243. This paper contains material presented 
at the probability seminar at the P.U.C. of Chile in September 2009. 

2. Notations and Main Results 

A Bouchaud trap model on Z with drift e, (X e (t); t £ [0, oo]) is a homogeneous Markov process 
with jump rates: 

c(x v)-( (l + e)r*-72iry = z + l ( ) 



where t = (t x ) x ^z are positive, i.i.d. under a measure P and satisfy 

lim u a P[r x >u] = l. (4) 
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For any topological space E, B{E) will stand for the a-algebra of Borclians of E. P x and E x 
will denote the probability and expectation conditioned on the environment r = (r x ) x£ z and with 
X e (0) — x. These probabilities are often referred as quenched probabilities. We define P x on 
Z N x R+ Z stating that for every A E B(Z N ) and B G B(R+ Z ), P X [A x B] := J B P x [C T }P(d T ), where 
C T := {x EZ N : (x,t) E Ax B}. 

P x is called the annealed probability. Note that X e is Markovian w.r.t. P x but non-Markovian 
w.r.t. F x . W is the expectation associated to P x . P° and E° will be simply denoted as P and E. 
Also P T and E T will stand for P° and E° respectively. These notations will be used with the same 
meaning for all the processes appearing in this paper. 

We have to make some definitions in order to state our main result: let B(t) be a standard one 
dimensional Brownian motion starting at zero and l(t, x) be a bi-continuous version of his local 
time. Given any locally finite measure fj, on R, denote 

Jr 

and its right continuous generalized inverse by 

VvO) := inf{s > : <^(s) > t}. 

The right continuous generalized inverse exists by definition, is increasing and, as its name 
indicates, it is a right continuous function. Then we define the speed measure change of B with 
speed measure X(/j,)(t) as 

X(ji)(t) := B(^(t)). (5) 
We also need to define speed measure changes of a drifted Brownian motion. Let C(t) := B(t)+t. 
We know that C(t) has a bi-continuous local time l(t,y). Given any locally finite measure fi in R 
we define 

<?V0) := / l(s,y)fj,(dy), 
Jr 

and its generalized right-continuous inverse by 

:= inf{s > : <^(s) > t}. 
Then wc define X(/j,)(t) (the speed measure change of C with speed measure /x) by 

X(ri(t):=C(Mt)). (6) 

By changing the starting point of our underlying Brownian motion B, we can change the starting 
point of X([i) and X(p). 

Let (xi,Vi) be an inhomogeneous Poisson point process on R x M+, independent of B with 
intensity measure av~ 1 ~ a dxdv. We define the random measure p as 

P-=^2viS Xz . (7) 
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The diffusion (Z(t);t £ [0,T]) defined as Z(s) := B(tp p (s)) is called the F.I.N diffusion. We 
also define the drifted F.I.N, diffusion Z(t) as Z(t) := C{xp p (t)). 

D[0,T] will denote the space of cadlag functions from [0,T] to K. (D[0, T],Mi), (D[0,T], Ji) 
and (D[0,T], U) will stand for £>[0,T] equipped with the Skorohod-Mi, Skorohod-Ji, and uniform 



Let V a be an a-stable subordinator started at zero. That is, V a is the increasing Levy process 
with Laplace transform E[exp(— AV^(t))] = exp(— t\ a ). Now we are in conditions to state the main 
result of this paper. 

Theorem 2.1. For all T > 0: 

(i) If a < a/(a + l) we have that {X^ n,a > {t)\ t £ [0,T]) converges in distribution to (V~ x (t);t £ 
(0,T)) in (D[0,T],U) where V~ l is the right continuous generalized inverse ofV a . 

(ii) If a = a/(a+l) we have that (X^ n - a \t); t £ [0,T]) converges in distribution to the drifted 
F.I.N, diffusion {Z(t);t£ [0,T]) on (D[0,T],U). 

(iii) Ifa>a/(a + l) we have that (X<- n ' a > (t);t £ [0;T]) converges in distribution to the F.I.N, 
diffusion {Z(t);t £ [0,t]) on {D[Q, T],U). 

We present heuristic arguments to understand the transition at a — jt^- First we analyze a 
sequence of discrete time random walks. Let (S e (i),i £ N) be a simple asymmetric random walk 
with drift e, S e (i) :— J2k=i b\i where is an i.i.d. sequence of random variables with: P[b£ = 

1] = l±i; p[6^ = -i] = We want to find the possible scaling limits of {S e{ - n \in);i £ [0,T]), 
depending on the speed of decay of e(n) to as n — > oo. 

We couple the sequence of walks S 16 ^) in the following way: Let (J7j)jgN be an i.i.d. sequence 
of uniformly distributed random variables taking values on [0,1]. We require that S e ^ takes 
his i-th step to the right {bf'^ = 1) if Ui > 1 ~^"- > and to the left otherwise. For each walk, 
we can decompose the steps into two groups: the first group is given by the steps i such that 
1 ~'^ n ' > < Ut < and the second group consists of the remaining steps. We can think that 

the first group of steps takes account of the drift effect and the second one takes account of the 
symmetric fluctuations of the walk. 

If the walk has given n steps, then the first group has about ne(n) steps, and the second group 
has fluctuations of order ^Jn. It is obvious that the drift effect will dominate the behavior if 



topology respectively. We refer to [14 for an account on these topologies. We define {X^ n ' a ^;t £ 
[0, T]), a rescaling of a walk with drift n~ a , by 
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i/n = o(e(n)). In this case we will have a ballistic (deterministic) process as a scaling limit. If 
e(n) = o(y / n) the fluctuations will dominate and we will have a Brownian motion as scaling limit. 
Finally the two behaviors will be of the same order if e(n) « i/n, and a Brownian motion with 
drift will be the scaling limit. 

The same reasoning can now be used to understand the change of behavior at a = a/ (a + 1) 
for the sequence of walks (X n (tn),t £ [0,T])„ £ n. In order to apply the precedent arguments we 
first have to estimate the number of steps that X n has given up to time Tn. To simplify we 
take T = 1. First, suppose that X n (n) is of order n u , where u is to be found. We know that 
after k steps, a walk with drift nr a is approximately on site kn~ a , so, it takes about n u+a steps 
to be on site n u . Thus, we can also deduce that at time n, X n has visited approximately n a 
times each site. As the distribution of satisfies (|3]), then the sum X)"=o Ti ^ s °f the same order 
that max{o<j< n u} n, and both are of order n u l a . We can estimate the time needed to arrive at n u 
as the depth of the deepest trap found (« n u / Q ) multiplied by the number of visits to that trap 
(sa n a ). This gives that n ss n = +a . But, we know, by definition, that X n arrives at the site n n l w 
approximately at time n. It follows that 1 = (u/a) + a, which yields u — (1 — a)a. This means 
that the number of steps that X n has given up to time n is of order n ( 1 - a )"+ a . 

Again, we can decompose the steps of X n into two groups. The first group accounts for the 
drift effect, and the second one accounts for the fluctuations. The first group will have approxi- 
mately n - a +[( 1 - a ) a + a } steps and the second group will give a contribution to the position of order 
n~ . Now it is easy to see that the ballistic behavior and the fluctuations will be of the same 

order i.f.f. [(1 — a)a + a]/2 = (1 — a)a or a — a/(l + a). 

3. The Supercritical Regime 

The proof for the constant drift case (a = 0) in [IS] is roughly as follows: first he prove that the 
sequence of rescaled first hitting times, (n -1 / Q inf{s > : X € (ns) > x} : x > 0), converges to an in- 
stable subordinator. Then, using that the right continuous generalized inverse of the process of first 
hitting times is the maximum of X e (t), he can deduce that (max{n~ 1 A e (n 1 / Q s) : s < t} : t > 0) 
converges to the inverse of an a-stable subordinator. Finally he shows that the walk and its 
maximum are close. 

For the proof of part (i) of theorem 12.11 we cannot follow the proof of [15] in a straightfor- 
ward way: suppose we show that a properly rescaled sequence of first hitting time processes 
(p a (n)H™(nx) : x £ R + ) (where p a {n) is an appropriate scaling) converges to an a-stable subor- 
dinator. Then, by inverting the processes, we get that the sequence (max{n _1 A™ (p a (n)~ 1 s) : 
s < t} : t £ R_|_) converges to the inverse of an a-stable subordinator. But we are searching a 
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limit for (m&x{d a (n)X n (tn) : t G K+}) (where d a (ri) is appropriate space scaling). That is, 
we want to obtain the limit of a sequence of rescaled walks where the drift decays as n~ a when 
the time is rescaled by n. But when we invert (p a (n)H"(nx) : x G M+), we obtain the sequence 
(max{n _1 X™ " (p a (n)~ 1 s) : s < t} : t > 0), which is a sequence of maximums of rescaled walks in 
which the drift decays as n~ a when the time is rescaled as p (n) _1 . 

To solve this, we will prove that the limit of (q a (n)H£, (nx) : x € R+) is an a-stable subordinator, 
where q a (n) is an appropriate scaling and b* sets an appropriate drift decay and depends on a. 
Inverting, we will obtain that (max{n _1 X" b (g a (n) _1 s) : s < t} : t > 0) converges to an a-stable 
subordinator. As we have said, we want the limit of a sequence of rescaled walks with a drift that 
decays as n~ a as the time parameter is rescaled by n. Hence, when the time parameter is rescaled 
as (7 (n) _1 , the drift should rescale as q a (n) a . Thus we need to choose b* so that n~ b = q a (n) a . 
But we know that q a (n) is the appropriate scaling for (H£,(nx) : x € R+). Hence, q a {ri) must be 
the order of magnitude of Hff,(n). That is q a (n) is of the order of the time that the walk X n 
needs to reach n. 

We now give a heuristic argument to find q a {n) and b* . When X n (t) has given k steps, 
it has an order kn~ b . So it takes about n b +1 steps to be on site n. We can think that the 
number of visits to each site x is evenly distributed. Then each site is visited about n b times 
before X n hits n. The time that the walks needs to reach n is of the order of the time spent 
in the largest trap. Thus we can estimate the total time spent by the walk as the depth of the 
deepest trap (which is of order n -1 /") multiplied by the number of visits to that trap. This gives 
a time of order n 1 / Q+h . What the previous arguments show is that X n b (t) arrives at n at time 
t w n 1 /" +b * (q a (n) w n 1 / a+f> *). But at that time we want to analyze a walk of drift {n 1 / a+b * )- a . 
That is, we need that a(l/a + b*) = b* . In this way we find that b* := a/[(l — a) a]. 

3.1. The embedded discrete time walk. For each natural n, the clock processes S n is defined 
as S n (0) :— 0. Furthermore S n (k) is the time of the fc-th jump of X n . S n is extended to 
all R + by setting S n (s) := S n ([s\). To each drifted walk I™' (t) we associate its corresponding 
embedded discrete time random walk (Y™ * :i£N) defined as Y" b := X n 6 (t) where t satisfies: 
S n (i) < t < S n (i + 1). 

Obviously Y™ b is a discrete time random walk with drift rT b . Wc can write 

fc-i 

i=0 

where (ej)j>o is an i.i.d. sequence of exponentially distributed random variables with mean 1. 
Define 

e = e(n) := n~ b 



BOUCHAUD WALKS WITH VARIABLE DRIFT 



9 



p = p(n) := (l + e(n))/2 

q = q(n) := (1 — e(n))/2 and 

v(n) := [c\og(n)n b J with c > 2. 

Let S(x, k) — 5(ar, fc, n) be the probability that y. e ^ hits x before fc starting at x + 1. Then 
we have that S(x, fc) = q+pE(x + 1, fc)S(x, A;) and that S(fc — 2, k) = q. These observations give a 
difference equation and an initial condition to compute S(x, k). Then we get that 

where r = r(n) := q(n)/p(n). Using that formula we can see that the probability that the walk 
ever hits x—1 starting at x is r. We now present a backtracking estimate. 

Lemma 3.1. Let A(n) := {mini<j<£ n ( n )(l^ n ^ — Y^ n ^) > —v(n)} where („(«) := minjfc > : 
Y* [n) = i}, then limn^ P[A(n)} = 1. 

Proof: We can write 

A c (n) = n \J{ min (Y« n) - ) < -„(«)} . 

Hence 



x=0 

But, in order to arrive from x to x — f(n), for each j = x — 1, . . . ,x — v(n), starting from 
j + 1 the random walk needs to hit j in a finite time. Hence, it takes v(n) realizations 

of independent events (strong Markov property) of probability r(n). In other words P[*4 c (n)] < 
nr ( n y( n ) = n {\ — 1+ ^ b « ) u ( n \ which can be bounded by n(l — ^pr) u ( n \ Replacing v{ri) we obtain 
n((l - ) cl °s(™). w e can see that (1 - ^pr)™ 6 -> e _1 when n -» oo. Now, for n big enough 

(l-^)" 6 * <e-i Then 

P[^ c (n)] < mi-5 c . 

But c > 2, so we get the result. 

□ 

Now we state the convergence result for the hitting time processes. 
Lemma 3.2. Let 

Then {H^(t)-t G [0,T]) convene* W eaA% to ((-^)-V<*V a (t);t G [0,T]) on (D[0,T],Mi), 
w/iere V Q (t) is an a-stable subordinates. 
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The proof of this lemma will be given in subsection 13.51 We present the proof of part (i) of 
Theorem 12.11 using lemma 13.21 and devote the rest of the section to the proof of lemma [ 



3.2. Proof of (i) of Theorem 12. li Let us denote 

X n (t) ^Ji^maxjr '*(snW a » +fc *);s€ [0,t]}. 

First we will prove convergence in distribution of X n to the (right continuous generalized) 
inverse of ( sin ^° a ) )~ 1 ^ a V a in the uniform topology. That is, we want to to prove convergence in 
distribution of the inverse of (ffW(t);t€ [0,T]) to the inverse of {{ si ^ a) Y 1,a V a {t)\t G [0,T]) in 
the uniform topology. Define 

C(T, S) n := {H^(S) > Tn^/ a ^ +b '}. 

Then, we have that, on C(T,S) n , the right continuous generalized inverse of (H^(s);s G 
[0,5]) is (X n (t);t G [0,T]). Let T > be fixed, by Lemma I3~2l we know that we can choose S 
big enough so that lim„_ i . 00 ¥[C(T, S) n ] is as close to 1 as we want. Let D^[0,T] be the subset 
of D[0,T] consisting of the increasing functions. By corollary 13.6.4 of [13], the inversion map 
from (D^[Q,T], Mi) to (D^[Q, T], U) is continuous at strictly increasing functions. Lemma (|3.2p 
gives convergence in distribution of (H^(t);t G [0,5]) to (( si ^" a) )~ 1/a V a (t);t G [0,5]) in the 
Skorohod M\ topology. We know that V a is a. s. strictly increasing, that is ({ sh ^ a ^ )~ 1 ^ a V a {t)\ t G 
[0,5]) G £>^[0,T] almost surely. So we can apply corollary 13.6.4 of [14] and deduce convergence 
in distribution of X n to the inverse of ( sin ^" a ) )~ 1 ^ a V a in the uniform topology. As we have 
said previously, the inverse of (H^(s);s G [0,5]) is (X n (t);t G [0,T]) in C(T,5)„. This proves 
convergence of the maximum of the walk. To deduce convergence of the walk itself it suffices to 
show that the walk is close enough to its maximum in the uniform topology. That is, to prove the 
theorem, it is enough to show that for all 7 > 0: 



sup In^X^ 

0<t<T 



{tn i - 1 /^ +b ')- X n {t)\ > 7 



Again, by Lemma l3~2l we know that F[Hg, (n log(ra)) > Tii' 1 ' 11 ^ 6 *] — > 1. Hence, we just have to 
prove that 



sup I n X r 

0<t<i?" (n log(n)) 



(t) - n- 1 m&x{X n (s); s G [0, t]}\ > 7 



-> 0. 



Which is to say, 



sup 

0<k<C n (ln\og(n)}; 



IV -Y k y >\ > 717 



where y c (") is the maximum of Y e ( n \ But, we can apply Lemma ETT1 to see that this is the case. 
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3.3. The environment. Here we give estimates concerning the environment. For each n € N 
define 

n l/a 

g( n ) '■— ~2~ - 

(Iog(n))i- 

Now, for each site x e N, we say that a; is an n-deep trap if > g(n). Otherwise we will say 
that x is an n-shallow trap. We now order the set of n-deep traps according to their position from 
left to right. Then call Si(n) the leftmost n-deep trap and in general call for j > 1, Sj(n) the j-th 
n-deep trap. The number of n-deep traps in [0,n] is denoted by 6(n). Let us now define 

£ l{ n) := {n^(n) (l - ^) < *„ < n v{ n) (l + ^ 
£ 2 (n) := {<5i A (mm (6, - Sj-i)) < p(n)}, 

l<7<0 n -l 

£3 fa) : = { max T ir < 5fa)}j an d 

— v(n)<x<0 

£(n) := £ifa) n £ 2 fa) n £ 3 fa) 
where p(n) :— n K k < 1 and <pfa) := P^ > fffa)]. 

Lemma 3.3. FFe have that lim^oc P{£(n)} = 1. 

Proof: 0fa) is binomial with parameters (n, <p(n)). £i is estimated using the Markov inequality. 
To control £2 it is enough to see that in 0, .., n there are 0(np(n)) pairs of points at a distance less 
than p(n). The estimate on £ 3 is trivial. 

□ 

3.4. Time control. In this subsection we prove results about the time spent by the walk on the 
traps. 

3.4.1. Shallow traps. Here we will show that the time that the walks spend in the shallow traps is 
negligible. 

Lemma 3.4. Letl(n) := j£i=o° T Y ^ e i l {r Y ,(n) < ^ "'togV)'"' }' Then 

P[Xfa)] -Jlasn^oo. (11) 

Proof: We have that P[Z(n) c ] = P[Z(n) c n£(n)] +o(l). Using the Markov inequality it suffices 
to show that 

/ n l/[(l-a)a]' 

l i>.'" ■ ■'" ^ ; = ° v log(n) 



E 



I] T ye( „)ea{r ( „ ) <9(„)}l {y «(")>_ l/( 



i=0 
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The number of visits of Y^ n ^ to x before time ( n (n) is 1 + G{x, n), where G(x, n) is a geometri- 
cally distributed random variable of parameter 1 — (q + pS(x, n)) (the parameter is the probability 
that, Y^ n \ starting at x, hits n before returning to x). Also 



Cn(n) 



Zl T Y^ e ' ll {r e(n )<9(n)} 1 



i=0 



(n) <<?(«)} i {y. £ (") >-v(n)} 



< ]T r x (l + E T [G(x,n)])l {Tsc<g{n)} . (12) 

x— — v{n) 



Using © we can deduce that (1 + E[G(x,n)]) < ^ ^"^ < cn b . So, averaging with respect 



to the environment in (|T2|) we get 



E 



C»(n) 



2j T Y i e( ' l ' e * 1 {T v e(») <9 W> 1 {V i e( ™ ) >-''(' 



i=0 



)} 



< Cn 1+b *E[r l {ro<g(rl)} ]. 



Also 



E[rol { . < 5 („) } ] < ^(l/2)^(n)P[r > (l/2)-'"+ 1 .g(n)]. 

Now, using (HJ) there exists a constant C such that the righthand side of the above inequality is 
bounded above by 

oo 

j=o 

Furthermore, since 1 — a > this expression is bounded above by Cg{n) l ~ a . This finishes the 
proof. 

□ 

3.4.2. Deep traps. Here we will estimate the time spent in deep traps. We define the occupation 
time for ifZas 

C»(") 

T x = T x (n) := ^ r^c^e; l {Y e( n)=x} . 

The walk visits x, G(x, n) + l times before Cn( n )> an d each visit lasts an exponentially distributed 
time. This allows us to control the Laplace transform of T x . For any pair of sequences of real 
numbers (a„)„ eN , (b n ) neN , a n ~ b n will mean that lim,,^,^ |^ = 1. 

Lemma 3.5. Let A > 0. Define X n :— 1/ ; ( ^_ a)Q ] . Then we have that 



E x [l - cxp(-\ n T x )\T x > g(n)} 



> 5(71)] -wa-" 



nsin(a7r) 

Proof: We must perform an auxiliary computation about the asymptotic behavior of the 
parameter 1 — (q + pS(x, n)) of G(x, n): 



(1 - (q +pE(x,n)))n b =p 



1 - r 
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2p(l + n~ b ) n ~ x 
_ (1 + n- b *)((l + n- b *) n - x - (1 - n~ b ') n ~ x ) 

= 2p 

(l + n-n(l-(l-TT^)"-*) 
which converges to 1. Thus we have showed that 



(1 - (q+pE(x,n)))n 



(13) 



We have 



E x [exp(-X n T x )]=E x 



exp 



G(x,n) 



where are i.i.d. end exponentially distributed with E(e,) = 1. Let A n := ^ a . Then 

E*[exp(-A n T x )] = 



1 + A„ 



n n b " (l-(g+pH(x,ra))) 

Using l[13|) we get that the above expression equals 

= = = 1 + o(n^ a ). 

1 + X„T X (1 + 0(1)) 1 + \ n T x 

Averaging with respect to the environment 

f°° 1 

E 2 '[l - cxp(-A„T 2; )l {r;c > g(n)} ] = / 1 - =— t 

J g(n) 1 + X n Z 



(dz) + o(n- 1 / a ) 



where the notation To(dz) denotes integration with respect the distribution of To. Integrating by 
parts J"^) 1 — z To(dz) we get that the above display equals 



X n z 

i + a; 



[to > z] 



g( n ) 



A, 



[t > z}dz + o(n~^ a ). 



!g(n) (1 + A„z) 2 

The first term is smaller than CA„(7(n) 1_Q = o(n ). To estimate the second term, note that 
for all r\ > we have 

(1 - ri)z- a < P[t > z] < (1 + r))z~ a 
for z large enough. Then we must compute f°? > - — x ~™ ^ z~ a dz. Changing variables with y = \ z 

to ° y J g(n) (1 + A„2) 2 to to " 1 + A„z 

we obtain 

K a ! . , , y- a (i-y) a d y . 

1 + A„s(n) 

But we know that this integral converges to T(a + l)r(a — 1) 



sin(7ra) ' 



□ 
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3.5. Proof of Lemma 13. 21 We will show the convergence of the finite dimensional Laplace trans- 
forms of the rescaled hitting times to the corresponding expression for an a-stable subordinator. 
This will prove finite dimensional convergence. 

Let = uo < • • • < uk < T and /3j, i — 1..K be positive numbers. We know that 



E 



K 



exp^-ft(( 



sin(7ra) 



r 1/a v a {ui) - ( 



7ra 



sin(7ra) 



K 



exp [J2 



sm(a7rj 



Uk-i. 



So, it only suffices to show that 

K 



E 



exp -fr(H (n) - H <"> (m-i)) 



K 



n— >oo 



exp > -— ; - u K - uk-i) 

\ Sin ttTT 



where is as in (fTOj) . We can decompose the trajectory of F £ (") up to Cn(L nu -ffJ) m t° three 

parts. The first one is the trajectory up to the time £ n (\_nuK-i ~ v (Tn)\), the second one is the 
trajectory between times ^ n (\jiUK-i ~ ^(Tn)^ and Q n {\nuK-i\ ), finally, the third part is the 
trajectory starting from time Cn(L nu Ar-iJ ) U P to time Cn(|_ nu -R"J )■ First we will show that the time 
spent in the second part of the trajectory is negligible. We have that F[m&x ye B lj(Tn) (x) > 9{Tn)] = 
o(l), which is to say that the probability of finding an n-deep trap in a ball of radius v{Tn) is 
small. Indeed Lemma 13.41 implies that there exists a constant C > such that 



t („)e. 



e,l 



(n) £ B ^(rn)(L 



"K-1"J)| 



< Cnd-)- (log(n))" 



Hence, the time that the walk spends in B v (r n )(\_UK-in\ ) is negligible. But in A(Tn) the 
walk never backtracks a distance larger than v(Tn), so, the time spent in the second part of 
the decomposition is negligible. The fact that in A(Tn) the walk never backtracks a distance 
larger than v{Tn) also implies that, conditional on A(Tn), the first and the third parts of the 
decomposition of the trajectory corresponds to independent walks in independent environments. 

So E[cxp(^^ 1 -Pi(HW(ui) - ffW^-i)))] can be expressed as 

K-l 



E 



exp -ft(ff (n) («.) - A" (n) K-i)) 



JgL««K-lJ 



exp 



-/M# (n) (^)-# Cn) K-i)) +o(l) 



where o(l) is taking account of the time spent in the second part of the decomposition of the 
trajectory and of A(Tn) c . 

The strong Markov property of y e( ") applied at the stopping time C, n (\nuK -ij ) and translational 
invariance of the environment give that H^(nui) — H^,(nui-i) is distributed as H£,(ns n (K)) 
where s n (K) — L"^ "J~^" jr - 1 "1 _ Iterating this procedure K — 2 times we reduce the problem to 
the computation of one-dimensional Laplace transforms. Hence, we have to prove that, for each 
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k < K 



,|exp 



(-0 k n-W°l- a H${na n (k)))] 



exp 



not 



sin(7ra) 



P% (ilk - u k - 



We have that ¥[S(Tn) H A(Tn)] -> 1, then we can write 

E[exp(-(3 k n-^- a H$,(ns n (k)))] = E[exp(-/3 fe n-( 1 / a )- a ffr.(^nW))l { £(T„)u^(T„) } ] + o(l). 

We know that the time spent in the shallow traps is negligible, so we only have to take into account 
the deep traps. We also know that on A(Tn), the walk does not backtrack more than u(Tn), and 
that, on £{Tn), the deep traps on [0,Tn] are well separated. Then we can write 



E[exp(-/3kn 



-{l/a)-a rrn 



Hi\(ns n (k)))}=E 



B(ns n (k)) 

J] EfMexp(-/3 fc n~ (1/a) - Q T 5! 
4=1 



o(l). 



Also, in £(Tn) we have upper and lower bounds for 9(Tn). Using the upper bound we see that 
the righthand side of the above equality is bounded above by 



E 



' n s n ( k ) tp ( n s „ ( k ) ) ( 1 - TZtfzkjkJ) ) 

JJ ^[expf-^n-^—TaJ] 



4=1 



o(l), 



Applying again the translational invariance of the environment and the strong Markov property 
we get that that the above display is equal to 

-(l/a)-arp s vyns n (k)tp(ns n (fc))(l- tog („^ n (fc» ) + 



E[E^[exp(-/3, 



which in turn can be expressed as 

E[exp(-/3 fe n-( 1 / tt )- Q r )|ro > g{ns n {k))] ns - (k)v(ns - {k)){1 ~ i^m) ) + (1). 
Using lemma p. 51) and the fact that s„(fc) — > — we obtain 



limsupE[exp(-^ fc n^( 1 / Q ^ a iJ™ (na„(fc)))] < 



exp 



a7r(3- a 



(u k - Uk-l) 



sin(a7r) 

The lower bound can be obtained in an analogous fashion. For the tightness, the arguments are 
the same as in Chapter 5 of [3] 



4. The Critical Case 

We want to show that for a — the sequence of walks (X^ n ' a \t); t € [0, oo]) converges in 
distribution to a drifted F.I.N, diffusion. We will mimic the arguments in [11]. But to treat the 
asymmetry of the model we will use a Brownian motion with drift instead of a Brownian motion. 
We use the existence of a bi-continuous version of the local time for a Brownian motion with drift. 
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4.1. The construction of the walks. Recall the definition of X(fi) given in display ((BJ. Let s 
be a real number and define 

\i := 2J ViS si . 
ie(z) 

Then X(fi) is a homogeneus Markov process with sL as its state space. The transition prob- 
abilities and jump rates of X(p) can be computed from the positions and weights of the atoms 
using the generator L of C(t) 

L "&:* . . ,14) 

The arguments we will give below are an adaptation of the reasoning used by Stone in [13j . For 
each i let r) si be the time of the first jump of X(n) started at si. By construction we will have that 
rjsi = Vil(a s ,0), where o~ s is the hitting time of (— s, s) by C(t). Using the strong Markov property 
for C(t) we can deduce that rj S i is exponentially distributed. It is easy to see that its mean is 
MjE[7(<7 SJ 0)]. Denote by Pt(x) the density at site x of the distribution of C(t) absorbed at {— s, s}. 
Using that l(a s ,0) := e^ 1 lim e ^ m(t € [0,ct s ] : C(t) £ [—e,e]) and applying Fubini's Theorem we 
find that E[l(a s , 0)] = eT 1 lim e ->o f° s P[C(i) G [-e, e]]dt. Then we find that 

E[/(tr„0)]= / pt(0)dt. 
Jo 

We also know that / °° pt(0)dt = /(0), where / is the Green function of (|14l) with Dirichlet 
conditions on {— s, s}. That is, / is the continuous function that satisfies 

+ fx = ~ So and /(s) = /(_s) = °- 

We know that the general solution to + |f = is g — C\ exp(— 2x) + Gi- This and the 

constraints on / give that 

ror 1 _ x exp(-2s) + l 

E IVsi] = v, L — — . 15) 

1 — exp(— 2s) 

For the computation of the respective transition probabilities we can use again the generator 
L. Let g : [— s, s] —> K be a continuous function such that + = an d g(~ s ) — 0, g(s) = 1. 
Using Ito's formula, we find that that g(C(t)) is a martingale. By the optional stopping theorem 
with the stopping time a s we find that the probability that the walk takes his first step to the 
right is 5(0). We can use the constraints on g to see that 

p[*(m)(t?«0 - s(i + 1)] = T . (16) 

1 + exp(2s) 

The proof of part (ii) of Theorem (|2.1|) will rely strongly on the following proposition. 

Proposition 4.1. Let (^ n ) n eN be a sequence of measures that converges vaguely to v, a mea- 
sure whose support is R. Then the corresponding processes (X(v n )(t),0 < t < T) converges to 
(X(u)(t),0 <t<T) in distribution in (D[0,T],U). 
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For the case where the underlying process is a Brownian motion, the proof of this fact can 
be found in [T3]. We will use the continuity properties for the local time I. For each fixed t, I 
is continuous and of compact support in x. Then, the vague convergence of v n to v implies the 
almost sure convergence of 4>i> n {t) to <j) v {t). As I is continuous in t, we obtain continuity of (j> Vn 
and of <f>w That, plus the fact that the <j) Vrl are non-decreasing implies that that <p Vn converges 
uniformly to <j> v . The function <j> v is almost surely strictly increasing, because the support of v is 
R. Now we can apply corollary 13.6.4 of [T3] to obtain that ip„ n converges uniformly to %j) v . That 
plus the continuity of the Brownian paths yields the lemma. 

4.2. The coupled walks. To prove part (ii) of Theorem 12.11 we will use Proposition 14.11 That 
is we want to show that each walk {X^ n ' a \t)]t G [0,oo]) can be expressed as a speed measure 
change of C(t), and then use convergence of the measures to get convergence of the processes. The 
problem is that we are dealing with a sequence of random measures, and the proposition deals only 
with deterministic measures. To overcome this obstacle we can construct a coupled sequence of 
random measures (p n )neN, such that (X(p n )(t);t € [0, oo]) is distributed as {X^ n ' a \t);t € [0, oo]) 
and that (p n )n<EN converges almost surely vaguely to p, where p is the random measure defined in 
(|7|) such that Z = X[p\. This section is devoted to the construction of the coupled measures. 

We recall that V a is an a-stable subordinator. To make the construction clearer, we will first 
suppose that To is equidistributed with the positive a-stable distribution V a (l). Let us consider 
the strictly increasing process (V p (t);t € R) given by V p (t) := p[0,t] if t > and V p (t) := —p[t,0) 
if t < 0. It is a known fact from the theory of Levy processes that V p {t) is a two sided a-stable 
subordinator. We now use this process to construct the coupled sequence of random measures 
(Pn)neN as 

i 

where s n := \ log -1 and 

rf := n 1 l^ a \V p (rT a ^ 1+a \i + 1)) - V p (n- a / il+a h)). (17) 

Observe that (r")igz is an i.i.d. sequence distributed like to, so that using (fTB"| and (fT5|) we 
see that X(p n ) is a walk with drift vT x l a taking values in s n Z. The latter means that X(p n ) is 
distributed like s„nTF? X^ n ' a \ The key observation here is that the scaling factor s n satisfies 

s n n a/{1+a) -^lasn^oo. (18) 

So, we just have to show that X{p n ) converges to X(p), because (fT8|) implies that if X(p n ) 
converges to X(p), so s n n a ^ 1+a ' > X(p n ) does. With (|T5)) in mind it is easy to prove that the 
sequence of measures (p n ) converges almost surely vaguely to p. Suppose that a < b are real 
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numbers and that V p is continuous at a and b, then 

Pn ((a,b]) = V p (n- a ^ 1+a ^a/s n \)~ V p (n- a ^ 1+a \[b/s n \ + 1)). 

But using (HE]) it is clear that [ a /s n \ a and \b/s n \ "4°° 6. Then the 

continuity of V" p at a and & implies that p n ((a, &]) n_ ij?° ^(a, 6], and we have proves the vague 
convergence of p n to p. 

Suppose now that To is not a positive a-stable random variable. Then, we can follow Section 
3 of [11] , There they construct constants c e and functions g e such that is distributed like t , 
where 

rt } := c-igdVM* + 1))-V p (d)). (19) 

Lemma 3.1 of [11] says that 

g e (y) -> y as 0. (20) 
As r satisfies (UJ) and using the construction of c e in Section 3 of [TT] , we can deduce that 

c e ~ e 1/Q (21) 

Define 

t™:^™^ (22) 

and again 

p n := ^ n -l/(l + a) r n 5sn .. 

i 

Then, by definition (TT5]), -^(pn) is a walk with drift n~ 1 l a taking values in s„Z. Using (j2"0"]l . 
(|2"TT) and (fT5)l we can see that P-a.s. p„ — > p vaguely. 

5. The Subcritical Regime 

We will prove that if a > a/(l + a), then {X^ n ^;t S [0, oo]) converges to a F.I.N, diffusion. 
We obtain the same scaling limit that was obtained in [11] for a symmetric B.T.M. Nevertheless, 
here we have to deal with walks which are not symmetric, in contrast with the situation of |llj . 
For this purpose we express each rescaled walk as a time scale change of a Brownian motion. The 
scale change is necessary to treat the asymmetry of the walk. Then we show that the scale change 
can be neglected. We now proceed to define a time scale change of a Brownian motion. Let /i be 
a locally finite discrete measure 

p,(dx) := WjS yi (dx), 

iGZ 

where (yi)iez is an ordered sequence of real numbers so that yi < yj i.i.f. i < j. 
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Let S : R — > K U {oo, —00} be a real valued, strictly increasing function, p will be the speed- 
measure and S the scaling function of the time scale change of Brownian motion. Define the scaled 
measure (So p)(dx) as 

(S op)(dx) := y^WjSs^^dx). 

i 

Let 

<t>(H,S)(t) := f l(t,y)(S o p){dy) 
and ip{fi, S)(s) be the right continuous generalized inverse of "0(a*> S). Then, as shown in |13] 

is a continuous time random walk with {t/j} as its state space. The mean of the exponentially 
distributed waiting time of X(S o p,) on yi is 

2m S {Vl+l )-S { y^) (23) 

and the transition probabilities to the right and to the left respectively are 

S(y i+ i) ~ S{y,) ^ S(yj) - S(y»_i) 



S(y i+ i) - S(yi-i) ' S(y i+1 ) - S(yi-i)' 
As in the previous section, we need to define a sequence of measures (v n )nez converging almost 
surely vaguely to p, and which can be used to express the sequence of rescaled walks X^ n ' a \ 
Let 

r n + 1 „ , 



—ft, — <, «,„•• ■ ; 



2r 



2n 
1 + / 

functions S n (which will converge to the identity mapping) given by 



where r" are defined in display (P22"j) , and r n := f — 1 T^ I -o ■ We will also use a sequence of scaling 



S»(i„«/(°+l)) :=V 

v ' n a/(a+l) 

We extend the domain of definition of 5™ to R by linear interpolation. Then, by ([23]) and 
we have that X(z/ n ,iS n ) is distributed like X^ n - a \ We will use the following theorem proved by 
Stone in [T3] . 

Proposition 5.1. Lei (f„)neN 6e a sequence of measures that converges vaguely to v. Then the 
corresponding processes (X(i/ n )(t), <t<T) converges to (X(^)(t),0 < t < T) in distribution in 
(D[Q,T],J X ) 



The proof of part (hi) of theorem 12. f I will rely in the following lemma. Let id denote the identity 
mapping on R, then we have that 
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Lemma 5.2. S n (n~ a '( 1+a ' \n a ^ a+1 ^> -J ) converges uniformly on compacts to id and v n to converges 
almost surely vaguely to p. 

Proof. The convergence of the scaling functions is easily seen to be true under the assumption 
a > a/(a + 1) because 

S n (n- a/{1+a) \n a/i - 1+a) x\) = 



n a/(a+l) 
j=0 



and 



n a/(l+a) — n a/(a+l) — n a/(l+a) 



Now we use the fact that 



1 + n~ a 



converges to 1, because a > a /(I + a). 

In a similar fashion it can be shown that the "correcting factors" in the definition of v n 

converge uniformly to 1 in any bounded interval. Hence, we can show the convergence of v n to p 
as in the previous section. 

□ 

Lemma 15.21 implies the vague convergence of (S n o v n ) to p. Then, by proposition 15.11 we 
can deduce that X(S n o v n ) converges to X(p). Let T > 0, by lemma I5T21 we have that S^ 1 
also converges uniformly to the identity. Thus, using the precedent observations, we get that 
(X(p,S)(t) : < t < T) converges to (X(p)(t)0 > t > T) in D[0,T] with the Skorohod J -I 
topology. We have proved that (X^ n ' a ^(t); t £ [0;T]) converges in distribution to the F.I.N, diffusion 
(Z(t);te [0,t]) on (D[0,T], Ji). 

Thus, it remains to prove that the convergence takes place also in the uniform topology. Using 
the fact that the support of p is K, we can show that 4>{p, id) is strictly increasing. The almost sure 
vague convergence of S o v n to p implies that, for all t > 0, (j>{v n , S n ){t) converges to <j>{p,id){t). 
As / is continuous in t, we obtain continuity of (j){y n , S n ) and of <p(p,id). That, plus the fact that 
the 4>(v n ,id) are non-decreasing implies that that <p(u n ,S n ) converges uniformly to <p(p,id). The 
function <p(p, id) is almost surely strictly increasing, because the support of p is R. Now we can 
apply corollary 13.6.4 of [2] to obtain that ip(v n ,S n ) converges uniformly to i/j(p,id). That, plus 
the continuity of the Brownian paths yields that X(S n ov n ) converges uniformly to X(p, id). Using 
that 5 1 " --1 converges to the identity, we finally get that X{y n , S n ) converges uniformly to X(p). 
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